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1. Introduction 

A quantum graph T is a graph (with multiple edges and loops allowed) in which each 
edge e is assigned a coordinate x (and hence a length) and the whole graph is equipped 
with a self-adjoint differential operator as Hamiltonian. For instance, if T is embedded 
nicely into a Riemannian manifold, one can use the arc length as a coordinate along an 
edge. In geometric language, T is a one- dimensional CW-complex with a Riemannian 
metric on each 1-cell and appropriate boundary conditions at each 0-cell to define a 
self-adjoint "Laplacian" . We refer to [8], [T2], [TJl 15] for further background on quantum 
graphs. 

Here we consider only graphs that are compact: the number of vertices V, the 
number of edges E, and the lengths of all edges are assumed to be finite. The number 
of edges attached to vertex v is called the degree of v and denoted d v . We also assume 
that every vertex has at least one edge attached, since isolated vertices are negligible in 
the quantum graph context. 

The simplest Hamiltonian for a quantum graph is the Laplacian with Kirchhoff 
boundary conditions, which acts as the negative second derivative along each edge, 

dXg ' 

with the functions in its domain required to be continuous at the vertices and to satisfy 
the Kirchhoff condition of no net flux at each vertex: 

e£E v 

Here E v is the set of edges incident on vertex v, and x e is the arc-length coordinate on 
e outgoing from v (in other words, the distance from v of the variable point on e). For 
more precise definitions see [121 US] an d section [3] below. 

In one of the earliest papers on quantum graphs [TS], J. -P. Roth calculated the 
trace of the heat kernel K for the operator just described. He found an exact formula 

oo „ 

yV An * = TrK = K(t,x,x)dx = Ki + K 2 + K 3 , (1) 

n=0 Jr 

where X n are the eigenvalues of H. Written in detail, it contains (before the integration) 
one term for every path in the graph T leading from the point x (not a vertex) to itself. 
These closed paths fall into three classes: The path of zero length yields the term 
K\ = L/yAni, the anticipated leading term in the Weyl series, where L is the total 
length of all edges of the graph. K 2 is the sum of the contributions of the periodic paths 
(where the initial and final direction of the path are the same), which are proportional to 
g-^c) 2 / 4 ^ where L(C) is the length of the path C. (Such terms do not contribute to the 
asymptotic expansion for t — > of Tr K in powers of t, but they determine oscillations 
in the distribution of the eigenvalues {A n }.) Finally, the contributions of paths that 
are closed but not periodic (i.e., the initial and final directions are opposite) sum to the 
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simple form 



K Z = \(V-E) 



(2) 



which constitutes the entire remainder of the Weyl series for the heat kernel trace. Note 
that K3 is independent of t and moreover is the only constant term in the formula (pQ). 

The expression in Equation <^ is interesting because it is a half-integer and depends 
on the topology of T only (e.g., it is independent of the edge lengths). Indeed, it is the 
Euler characteristic of T regarded as a 1-complex. These features are reminiscent of 
index theorems in geometric analysis and the calculation of indexes from the constant 
terms in heat-kernel expansions [6l [7J. The original goal of this paper was to give an 
index interpretation of ([2j); in fact, we also generalize it to graphs with other boundary 
conditions and compute indexes of quite general quantum graph operators by another, 
very simple method. 

The contents of the paper are as follows: Section previews the appearance of indexes 
in the heat kernel asymptotics for the case of an interval. In section [3] we introduce 
necessary notions and auxiliary results concerning quantum graphs. The section also 
contains a general formula for indexes of differential operators on quantum graphs. 
Section H] contains the main results concerning the relations between the constant terms 
in the asymptotic expansion of the heat kernel and indexes of suitable operators on 
the graph. Relations to the secular determinant are discussed in section [5j The final 
section [6] contains some remarks and conclusions. 

2. The interval 

In [7J, Section 1.5] P. Gilkey treats the Laplacian on an interval with Dirichlet and with 
Neumann boundary conditions as the prototype of the index theorem for the de Rham 
complex on a manifold with boundary. The index theorem for a quantum graph with 
Kirchhoff boundary conditions is a different generalization of this elementary example, 
so we shall review the latter. 

Let Hp and be the operator — d 2 /dx 2 on the interval (0, L) with Dirichlet and 
Neumann boundary conditions, respectively. The eigenfunctions of with eigenvalue 
are the constant functions, so the kernel of Hm has dimension 1. In contrast, Hp has 
trivial kernel. On the other hand [7J Subsection 3.1.3], the heat traces of these operators 
are 



where the negative sign applies to the Dirichlet case and the exponentially small terms 
analogous to K2 in (CQ) have been omitted. Therefore, 



To identify (j4]) as an index theorem we must factor H^ ^ into first-order operators. 
Let A be the operator d/dx acting on the domain H l (0, L), which is the Sobolev space 




(3) 



dim ker Hn — dim ker Hp = 1 = Tr Kjq — Tr Kd . 



(4) 
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containing functions on (0, L) that, together with their first distributional derivatives, 
are square-integrable. Standard integrations by parts show that the adjoint operator 
A* is — d/dx with domain Hq(0, L) (containing functions from i? 1 (0, L) that satisfy the 
Dirichlet conditions /(0) = f(L) = 0) and that A** = A. One now forms two second 
order operators 

H N = A* A, H D = AA*, (5) 

where in the first case the domain consists of twice differentiable (H 2 (0,L)) functions 
satisfying the Neumann conditions, /'(0) = f'(L) = 0, so that Af G domA* and the 
composition is defined; in the second case, similarly, the domain consists of functions 
from H 2 (0,L) satisfying the Dirichlet conditions. 

Because (f,A*Af) = (Af,Af) = \\Af\\ 2 , the kernel of A is the same as that of 
if at . Similarly, ker A* = ker Hp . Therefore, (j3J) can be restated as the index formula 

index A = Tr K N - Tr K D = 1. (6) 

Remark 1. By regarding the elements of domA* as 1-forms rather than scalar 
functions, one identifies A and A* with the exterior derivative operator d:A°(0, L) — > 
A\0,L) and its adjoint 5:A\0,L) -> A°(0,L), where A°(0,L) and A^O.L) are the 
spaces of L 2 functions f(x) and 1-forms g(x) dx. One can combine these operators into 
a single operator from A°(0, L) © A 1 (0, L) into itself, 

where, in the version on the left, d annihilates the 1-forms and S annihilates functions. 
Then the Hodge Laplacian is 

(d + 5) 2 = 5d + d5 = 

It is this formulation that generalizes to higher-dimensional manifolds, with Hn acting 
on forms of even degree and Hp on forms of odd degree (or vice versa) [HI [7j. 

In the following sections we will extend this analysis to more general quantum 
graphs. In particular, a central task is to identify the analogues of the operators A 
and A*. 

3. Quantum graphs 

3.1. Vertex conditions 

As we have mentioned in section appropriate vertex conditions are needed in order 
to turn the (negative) second derivative along the edges of a quantum graph into a 
self-adjoint operator in L 2 (T). All such choices of boundary conditions at vertices were 
catalogued in [12] (after prior discussion in [1]) and reformulated in [9l [15]. It will be 
convenient for us to follow the formulation from [15J. 




Index theorems for quantum graphs 



5 



Let v be a vertex and f(x) a function on T. We denote by 



/ fi{v) \ 



F(v) 



\ /*» J 

the vector of values of the function / at the vertex v, attained along d v edges incident 
to v. In particular, if / were continuous, all these values would be equal. Analogously, 



F'{v) 



is the vector of derivatives at v of / along these edges, where the derivatives are taken 
in the directions outgoing from the vertex v. 

It is clear that vertex conditions for the second-derivative operator can involve only 
the values of the function and of its derivatives along edges. If these conditions do 
not mix the values attained at different vertices, they are called local. (On an interval, 
for instance, Dirichlet, Neumann, and Robin conditions are local, but the periodicity 
condition is nonlocal.) As we will see later, there is actually not much difference between 
local and nonlocal vertex conditions on a quantum graph. (For instance, the periodicity 
condition becomes local if the interval is replaced by a loop attached to a single vertex.) 

Theorem 2. [15] All self-adjoint realizations H of the negative second derivative on 
T with local vertex boundary conditions can be described as follows: For every vertex 
v, of degree d v , there are two orthogonal (and mutually orthogonal) projectors P v , Q v 
operating in C dv and an invertible self-adjoint operator A v operating in the subspace 
(1 — P v — Q v )<C dv . (Either P v , Q v , or C v = 1 — P v — Q v might be zero.) The functions 
f in the operator domain are those members of the Sobolev space ® e H 2 (e) that satisfy 
at each vertex v boundary conditions consisting of the "Dirichlet part" 

P v F(v) = 0, (7) 

the "Neumann part" 

0, (8) 



QvF'(v) 
and the "Robin part" 

CF'M 



KC v F(v). 



The quadratic form of H is 



h[f,f]=J2[ J- dx + Y,{KC v F,C v F) 



(9) 



(10) 



with the domain that consists of the functions f(x) that belong to the Sobolev space 
if 1 (e) on each edge e and satisfy §7ty at each vertex. 
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Remark 3. This theorem was formulated a little bit differently in [15] . For one thing, 
was called — L v there. More importantly, there the two projectors Q v and C v where 
lumped into a single one, and thus the condition of invertibility of the operator 
disappeared. The equivalent formulation provided here distinguishes between Robin 
and pure Neumann conditions, as is often convenient. 

Remark 4. The three parts (ED)-© of the vertex conditions can be combined into a 
single condition 

A v F{v) + B v F'{v) = 0, (11) 

where the d v x d v matrices A v and B v are 

A v — P v — A V C V , B v = Q v + C v . (12) 

The conditions were introduced in [12] in the form ( TTTT) (which by itself does not suffice 
to define the matrices A v and B v uniquely, however). 

The most popular vertex conditions are the Kirchhoff ones (also called Neumann 
or natural), which reduce at vertices of degree 1 to Neumann conditions: 

Definition 5. The Kirchhoff boundary conditions are defined by the continuity condition 

fi(v) = Mv) = - = f d .{v) = f(v) (13a) 

as Dirichlet part and the zero flux condition 

£i» = o (135) 

e=l 

as Neumann part, with no Robin part. 

In other words, (1 — P v )C dv is in this case one-dimensional and consists of the 
vectors with equal coordinates. 

Another type of conditions that arises in our work is dual to the Kirchhoff type, in 
the sense that the roles of the values and derivatives of the function / at each vertex 
are switched. (At vertices of degree 1 these "anti-Kirchhoff" conditions reduce to the 
Dirichlet ones.) 

Definition 6. The anti-Kirchhoff boundary conditions are 

d v 

X>M = (14a) 

e=l 

as Dirichlet part and 

m = m = ... = f dv {v) = f{v) (146) 

as Neumann part, with no Robin part. 
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3.2. Bonds vs edges 

In what follows, we will need to use directed edges (which we will call bonds) rather 
than the undirected ones as before. Thus, each edge results in two directed bonds (with 
opposite directions), which are denoted by Greek letters. We denote by a the bond a 
with its direction reversed. 

Recall that loops (tadpoles) can always be removed from a quantum graph by 
inserting extra Kirchhoff vertices of degree 2. Adding such a vertex does not change the 
heat trace or the Euler characteristic, nor either side of any of the index formulas in this 
paper. Therefore, one may assume that the two ends of a bond are distinct vertices. 

It is not necessary to pick either of the two directions of an edge as the canonical 
one. The language of differential forms makes it possible to give global meaning to 
the differential of a function on T without committing to any particular coordinate, 
x e , on each edge. In discussing the behavior of functions (and their derivatives) in the 
neighborhood of any one vertex, therefore, we remain free to use the most convenient 
coordinate on each edge, namely, the outgoing arc length parameter. 

3.3. Scattering matrices 

In this subsection we introduce, following [12l E] , the scattering matrices and some of 
their properties that we will need in the rest of the text. 

Let if be a self-adjoint realization of the negative second derivative — d 2 /dx 2 on 
a finite quantum graph T (i.e., one of the self-adjoint vertex conditions described in 
Theorem [2] is imposed) . 

Let us consider a vertex v and the set E v of all edges e incident to it. (Such a 
configuration is called a star] see Fig. CD) For any edge eo G E v and any real k, we 



choose as in section [T] the coordinate x increasing away from the vertex and consider 
the unique solution f(x) on the star E v of the following scattering problem at v: 



In other words, o"e e is the reflection coefficient along the bond eo , and ae i is the 
transmission coefficient from the edge eo to e. Notice that the coefficients ai^l in general 
depend on k. 




Figure 1. A star. 




V ) 



(15) 
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Definition 7. The unitary d v x d v matrix a^(k) with the entries ai% 2 (k) for ej G E v 
is the edge scattering matrix at the vertex v. 

Notice that in defining the direction chosen along each edge depends on the 
vertex considered. That is why it becomes necessary to deal with directed bonds when a 
scattering matrix S for the whole graph is defined. However, as explained in section |3~2| 
the ambiguity in x e does not create any inconsistency in the notation. Another remark 
is that this matrix clearly depends upon what type of vertex conditions are imposed, and 
not every matrix function a^ v \k) can necessarily be realized by one of the second-order 
differential Hamiltonians studied here. 

It is straightforward to derive the formula [12] 

cr (v) (k) = —(A v + ikB v )-\A v - ikB v ), (16) 

which in particular confirms that the matrix is unitary and shows that its /^-dependence 
is tightly constrained. From (jl2j) we get an alternative representation of a: 

a^(k) = P V -Q V + (A v - ik)-\K + ik)C v . (17) 

The following result will be important for what follows. A part of it was proved by 
Kostrykin and Schrader ( [TlT Proposition 2.4]), [T2J, Corollary 2.3], [T3J Theorem 1]). 

Theorem 8. The following conditions are equivalent: 

(i) For each vertex v, <r^ v ' is independent of k. 

(ii) For each vertex v, there is a value k ^ such that (er^) 2 = 1. 
(Hi) For each vertex v, (a^) 2 = 1 for all k. 

(iv) For each vertex v, has the form 1 — 2Q V for some orthogonal projection Q v . 

(v) There is no Robin part in the vertex conditions: C v = for each vertex v. 

(vi) The vertex conditions are scale-invariant (i.e., if a function f(x) on neighborhood 
of v in the star E v satisfies the vertex conditions at v, then after rescaling to f{rx), 
it still satisfies the conditions) . 

(vii) The Hamiltonian H = — with the given vertex conditions can be factored as 
H = A* A, where A = with appropriate vertex conditions, and A* is its adjoint 
operator. 

Proof. Equivalence of statements (i) through (iv) was proved by Kostrykin and Schrader; 
it also follows rather easily from ([T7|) . 
Computing (o^) 2 using ( 1T71) . we get 

(° {v) ) 2 = p v + Q*> + (t~^) 2 Cv ■ (18) 

If C v 7^ 0, then in order to get a 2 = 1, we need that ^ ^-ifc ) = 1- A straightforward 
calculation shows that this is impossible for an invertible operator A and a non-zero k. 
This proves equivalence of (iii) and (v). 
Equivalence of (v) and (vi) is trivial. 
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The implication (vii) =>- (v) can be established as follows. The vertex conditions 
for A can involve only the values of the function, not its derivatives. Thus, at any 
vertex v they can be written as P v F(v) = for some orthogonal projector P v . Then 
a simple and well known (e.g., [12]) calculation, which boils down to an integration by 
parts, shows that for the adjoint operator — -r- , the vertex conditions are given by the 
complementary projector Q v — 1 — P v . The equality H = A* A now implies that the 
vertex conditions for H have P v as the Dirichlet and Q v as the Neumann part, with no 
Robin part being present. This argument can easily be reversed to show the converse 
implication, (v) =^> (vii). □ 

Corollary 9. The scattering matrices a for the Hamiltonian — -r4 with Kirchhoff or 
anti-Kirchhoff boundary conditions satisfy the equivalent conditions of Theorem 

We now introduce the global scattering matrix S a @ , entries of which are indexed by 
the directed bonds a and f3. 

Definition 10. The 2E x 2E (global) scattering matrix S is defined as follows: 

{Ual if ol terminates at v and 8 starts at v, 
Pa J (19) 
otherwise, 

where a and (3 are (directed) bonds in V. 

Proposition 11. Let H satisfy the equivalent conditions of Theorem^ so that it factors 
as H = A* A as in (vii) of the theorem. Let also H' = AA* . We also denote by a' and 
S' the scattering matrices for H' . Then 

3 a /3 — —{3 )a/3 

for any vertex v and any bonds a, 8. 

Proof. Indeed, it is clear that for H' the projectors P v and Q v exchange their places, 
while C v — 0. Then formulas (fTTj) and f|T9l) prove the statement. □ 



3.4- Indexes of quantum graph operators 

As it happens, one can establish a simple formula for the index of any (elliptic) 
differential operator on a compact quantum graph T, which in particular implies the 
index formulas for the exterior-derivative operators A introduced previously. 

First we need to review the basic notions concerning the Fredholm property and 
the index (see, e.g., [I0])- Recall that the codimension of a (closed) subspace E C H 
is defined as the dimension of the quotient space H/E, or, equivalently (in a Hilbert 
space), the dimension of an orthogonal complement of E. 

Definition 12. A bounded operator T: H\ — > H2 between two Hilbert (or Banach) spaces 
is said to be Fredholm, if it has a closed range and the dimension of its kernel kerT 
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and the codimension of its range ran T are finite. The index of a Fredholm operator T 
is defined as 

index T = dimkerT — codimranT. 

Proposition 13. If operator T is Fredholm and operator K is compact (in particular, 
of finite rank), then T + K is also Fredholm and index(T + K) — index T. 

To formulate the main theorem for a differential operator of arbitrary order it is 
convenient to choose an orientation for each edge, so that the arc length parameter x e 
is unambiguous. 

Theorem 14. Consider the operator on T defined by the differential expression of 
order m 

with Cq(x) continuous on each closed edge (but not necessarily on the whole graph) and 
never equal to and all other Cj measurable and bounded. Let T\ be the restriction of 
T as an operator from © e iJ m (e) into L 2 (T) to a subspace of codimension p (e.g., by 
imposing p vertex conditions sustainable by H m , i.e., involving derivatives up to the 
order m — 1). Then 

(i) The operator so defined is Fredholm. 

(ii) index 7\ = mE -p. (22) 

Proof. Consider T as the naturally defined (and obviously bounded) operator from 
© e if m (e) into L 2 (T). All terms in T that involve derivatives of order less than m are 
compact operators and thus do not influence the Fredholm property or the index. Thus, 
we can assume that T = c (x e )£^;. This is now the composition of acting from 
@ e H m (e) to L 2 {Y) with the invertible operator of multiplication by c (x) in L 2 (Y). Thus, 
everything reduces to the mth derivative alone. It is easy to show that it is a surjective 
operator from H m (e) onto the whole L 2 (e) (and thus from @ e H m (e) onto L 2 (T)). On 
each edge, it has the m-dimensional kernel consisting of polynomials of degree less than 
m. Thus, T is Fredholm and 

index T = J^(m - 0) = mE. 

e 

Let us now notice that by definition, 7\ is the restriction of T onto a subspace M 
of codimension p. Consider any (p-dimensional) complement iV to M in © e iJ m (e) and 
the extension T of T\ from M to the whole © e if m (e) that acts as the zero operator on 
N. Then the difference T — T vanishes on M and thus is a finite-dimensional operator. 
Hence, T is Fredholm of the same index mE as T. 

On the other hand, it is clear that the ranges of T and Ti are the same and the 
kernel of T is p dimensions larger than the kernel of T. Hence, T is Fredholm and 
index T = index T — p = mE — p. □ 
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This implies in particular 

Corollary 15. Let A be the exterior derivative d acting as a bounded operator from 
a subspace M of codimension p in @ e H 1 (e) into L 2 (T). This operator is Fredholm of 
index E — p. Therefore, 

indexA = E -p = E -^2<limP v , (23) 

where P v are the orthogonal projectors describing the vertex conditions for A. Thus, in 
particular, 

(i) Without any vertex conditions, one has index A = E. 
(ii) With continuity conditions $13a\) at all vertices, one has 

index A = E - - 1) = E - (2E - V) = V - E. 

V 

(Hi) With the condition |7Ja|] that the sum of values at each vertex is equal to zero, one 
has 

index A = E - ^1 = E - V. 

V 

Notice that ^2 v dimP v arising in this corollary is just the number of vertex 
conditions defining H that contain only the values of the function and no derivatives 
(Dirichlet part of the conditions). 



4. Heat kernel and index 



In this section we address the relation between the heat trace asymptotics and the index 
on quantum graphs. Most of the considerations are independent of Theorem [TU 

We will assume from now on that V is an arbitrary quantum graph and the 
Hamiltonian H satisfies the conditions of the Theorem [HI so that it factors as H = A* A, 
where A = d/dx with some vertex conditions on the values of functions, such conditions 
corresponding at any vertex v to an orthogonal projector P v in C Then, as before, 
we denote by H' the operator AA* with the vertex conditions given by the orthogonal 
projector Q v = 1 — P v . Let also K and K' be the corresponding heat kernels. The 
following proposition is standard. 

Proposition 16. In the situation just described, 

index A = Tr K - Tr K'. (24) 

Proof. The non-zero eigenvalues of H = A* A and H' = AA* are the same, including 
their multiplicity. The only exception is that the dimensions of the eigenspaces for the 
eigenvalue are different. Thus, the difference of the heat kernel traces is guaranteed 
to be independent of t. At large t this difference reduces to the difference of the 
nullities (i.e., dimensions of the kernels), and at small t it reduces to the difference 
of the constant terms in the heat-kernel expansions. See, for instance, [6] for a more 



Index theorems for quantum graphs 



12 




Figure 2. A bounce path aa^a-ia\a with n = 3. 



detailed exposition. Now, since the nullity of H is clearly equal to that of A, and that of 
H' coincides with that of A*, one concludes that Tr K— Tr K' = dim ker if —dim ker H' = 
dim ker A — dim ker A* = index A. □ 

We now need to establish a formula for the constant term of a heat trace: 

Theorem 17. [HI [1] Let T be a finite quantum graph with scale-invariant vertex 
conditions defining the Laplacian. (Thus the bond-to-bond scattering matrix S a p is 
independent of the frequency k). Let K(t,x,y) be the corresponding heat kernel on T. 
Then the constant term in the asymptotic expansion at t — > of the heat trace 

oo „ 

^V A "* = / K(t,x,x)dx (25) 



n=0 

is 



i5> 



,^ .... (26) 

a 

(See Definition IW\ for Sap.) 

Sketch of proof. The theorem is proved in [1] for a different kernel, but as stressed in 
[19] the same argument applies to a whole class of kernels, including the heat kernel. 
(See also [TBI ITT] .) Starting from the standard one-dimensional heat kernel on the real 
line, 

K (t, x, 0) = (4vrt)- 1 / 2 e- a;2/4 *, (27) 

by an extension of the method of images one constructs the heat kernel on the graph 
as a sum over all paths from y to x, which then needs to be restricted to the diagonal 
y = x. The heat trace is formed then by integrating over x. As in [18], the contributions 
of the periodic paths, i.e., the ones that return to the point x with the same direction 
as at the start, are proportional to Gaussian terms e _i ^ //4t and cannot contribute to 
the t-independent term of the heat-kernel expansion. The path of zero length gives 
the leading Weyl term, proportional to t~ x l 2 . The contribution of the other class of 
paths, which are closed but not periodic ("bounce" paths), apa, where p is a cycle in V 
(see Fig. [21 where the triangle represents the cycle p and the point x is located on the 
bond a), can be reduced to the following sum: 

a 



I 

n peP n 



2 X) Z) Z) Aa ^ / x) dx - (28) 
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(The condition of ^-independence of the scattering matrix is used here.) We denote here 
by P n the set of cycles of period (number of edges traversed, including multiplicity) n. 
E.g., in Fig. El the cycle p = aic^a^ has period 3. We also use here notations L a for the 
length of the bond a and / p for the metric length of the cycle p (e.g., l p = L ai +L a2 +L a3 
in Fig. E])- The shorthand notation A ap a is used for the product of scattering amplitudes 
along the path 

A _ — q q ...£?_ 

One can now make a sequence of reductions [T^J |TJ: The unipotency property, a 2 = 1, 
of the scattering matrix (see Theorem [8]) leads to 

a 

This allows a massive inductive reduction of the sum fl28|) . in the course of which one 
must also combine the a-dependent factors J Z P a Ko(t,x) dx. One eventually arrives 
at the following representation of (1281) : 



1 f°° 

-J2 S ™ K (t,x)dx. (29) 

From (1291) and (|27|) one obtains (1261) as the total contribution of all these "bounce" 
paths. This finishes the proof of Theorem [171 □ 

Example. It is well known (e.g., [H]) that at a Kirchhoff vertex of degree d v the 
scattering matrix is 

4 ] - T - <W • (3°) 

(l v 

For an entirely Kirchhoff graph, therefore, we have 



4 ^ V d v J 2 2 

a v=l x ' 

since every edge is incident on two vertices. This reproduces Roth's formula (EJ). By 
a similar calculation, or by appealing to Proposition [TTJ one sees that in the case of 
a graph all of whose vertices are of the anti-Kirchhoff type, the resulting term is the 
negative of fl3T|) . 

Proposition 18. In the context of Proposition [TR Tr K — Trf^' is equal to twice the 
constant term in the small-t expansion ofTrK. 

Proof. Proposition [TT] and Theorem [T7] imply that the constant term in Tr K is the 
negative of that in Tr K' . The proof of Proposition [16] shows that only these constant 
terms survive when the traces are subtracted. □ 
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Corollary 19. Under the conditions of Theorem^T^ and Proposition 1Tb} the following 
alternative representations hold for the index of A: 

mdexA= - S a s = Tr K - Tr K' = E - ^ dimP,, — E — p, (32) 

a v 

where H = A* A is the factorization of the Hamiltonian in accordance with (vii) in 
Theorem^ K and K' are the heat kernels of H and A* A, E is the number of undirected 
edges in V , P v is the projector onto the Dirichlet part of the vertex conditions, andp is the 
total number of vertex conditions not involving derivatives. Furthermore, Tr K — Tr K' 
can be read off from the asymptotics of a single heat kernel by virtue of Proposition [TR 

Proof. The last two equalities are quoted from ( l23l) for completeness. The rest of the 
corollary summarizes the results of this subsection. □ 



4-1. The Euler characteristic 

We now look at the special situation of Kirchhoff conditions to see the implications of 
Corollaries [19] and [15] there. 

It is convenient to consider first-order operators on quantum graphs as defined in 
terms of differential forms, rather than functions, introducing thus an analogue of the 
de Rham complex. Therefore, we henceforth identify A with d: A°(r) — > A : (r) and A* 
with 5: A X (T) — > A°(r) (on appropriately restricted domains). 

Theorem 20. Let Hk and Ha be the Kirchhoff and anti-Kirchhoff Laplacians on a 
compact quantum graph T (acting in A°(r) and A x (r) respectively), Kk and Ka be the 
corresponding heat kernels, and d and 5 = d* be the external derivative operators ( with 
the domains defined by the continuity conditions for d and the sum of values equal to 
zero at each vertex for 5 ), so that 

H K = 6d and H A = d8. (33) 

Then 

index d = Tr K K - Tr K A = V - E, (34) 
the Euler characteristic ofT. 

Proof. Either the first formula (with (13 ip ) or the last formula in (1321) can be used to 
calculate the index. □ 

Corollary 21. Let C be the number of connected components ofT. Then 

dim ker d = dim ker Hk = C (35) 

and 

dim ker 5 = dim ker Ha = E — V + C. (36) 
In particular, in the connected case 

dimker5 = P-1/ + l = r, (37) 
the rank of the fundamental group ofT. 
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Proof. Equality (1351) is immediate, since the zero modes of are constant on each 
component. Now the index theorem ([341 yields fl36|) . □ 



Remark 22. Equation f[36[) counts the locally constant differential 1-forms satisfying 
(1 14 aft . It is for a tree graph (where 1 1 14 aft must be violated at the leaves) and 1 for a 
cycle. In general, it counts the independent cycles in the graph. 

5. Relation to the secular determinant 



Kottos and Smilansky [14] derived their trace formula for the density of states of a 
(compact) quantum graph from a certain secular equation, 

f(k) = det [[/(&) - 1] = 0. (38) 

In terms of frequency the density of states is 

oo 

p (k) =J2 5 ( k ~ \ k n\) = N S(k) + Pl (k), (39) 

n=0 

where iVo is the true spectral multiplicity of as an eigenvalue, and pi(k) is the 

contribution of the positive eigenvalues, A n = fc n 2 , appearing with multiplicity. (We 

assume there are no negative eigenvalues. According to [HI [TFl El [16], for strictly 

positive eigenvalues the spectral multiplicity is equal to the multiplicity as a root of 

( [381 . but the situation for k = is quite different.) The secular function / defines a 

distribution p(k) on the entire real line by 

r„ i 1 f°° ff'(k-ie) f(k + ie)\ ,,,, „ , , 

p[4>] = lim — / ^- — - -rn — rf 4>{k) dk 40 
40 Am \ f(k - le) f(k + ie)J^ J 

oo 

= 2 E <t>(k n ), (41) 

n=— oo 

the sum being over all zeros of / (with multiplicity), including the possible one at k = 0, 
whose algebraic multiplicity as a root of (|38[) we shall denote iV. If /c is a nonzero root 
of /, then so is —k. The positive eigenvalues are the squares of the nonzero roots of /; 
thus every positive eigenvalue appears in the sum ([411 twice (times its multiplicity), and 
hence ( 1391) is correctly reproduced on the positive axis by p. The spectral multiplicity 
Nq , however, is generally equal neither to iV nor to ~iV. Furthermore, N itself has 
proved difficult to calculate reliably [HI [13, [HI EE] • 

From the foregoing definitions and discussion it follows that 

(iVo - ~N)8(k) + p(k) for k > 0, 

for k < 

= /i(A0[(A r o-|iV)<H£0+/W], (42) 
where h is the unit step function and (because the product hS would otherwise be 
ambiguous) we stipulate that h6[<f>] = 0(0) (not |0(O)). In terms of the density of 
states, (p]) can be rewritten as 



p(k) 



e' k l p(k) dk = Tr K = K x + K 2 + AT ; 



3 • 
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The derivation of the trace formula [T4"l [T7J [SJ [TB] makes clear that p yields precisely the 
leading Weyl term and the oscillatory (periodic-orbit) terms in the density of states. It 
follows (cf. discussion following (DO)) that the contribution of p to the heat kernel is the 
terms K\ and K 2 . The remaining term in fH2|) . proportional to S(k), must therefore be 
responsible for precisely the constant term K3 in the heat kernel. We reemphasize that 
the coefficient of those terms is not N , the coefficient of 5(h) in the spectral density p. 
That is not a paradox: The distributionally convergent periodic-orbit sum p contains 
another 5(k) contribution to p that restores consistency with ( l39i) . 

On the other hand, we now know from previous sections that K 3 , for a scale- 
invariant Laplacian, is equal to half the index of the associated first-order operator, A. 
The index, in turn, equals E — p, where E is the number of (undirected) edges and p is 
the number of Dirichlet conditions. Furthermore, we have index A = Nq — Nq, where 
Nq is the nullity of A* and hence of the operator dual to H. Therefore, we immediately 
get two interesting identities: 

Corollary 23. Let Nq and N be the spectral and algebraic multiplicities of k = for 
a scale-invariant graph Laplacian, A* A, and let Nq be the spectral multiplicity for the 
dual Laplacian, AA* . Then 

N = 2N - index A = 2N - E + p (43) 

and 

N = N + N*. (44) 

Example 1: For Hk , the Laplacian of a connected Kirchhoff graph, one has Nq — 1 
and index A = index d — V — E. Therefore, N = 2 — V + E, in agreement with [T7J 
corrigendum] and [IB] . 

Example 2: For the pure Neumann Laplacian of section 13.41 one has N = E 
and p = 0, so iV = E. This is correct, because appears as a root of / once for each 
disconnected Neumann edge. 

Kurasov [16J gives a convincing direct calculation of N for the Kirchhoff case. On 
that basis he deduces that the Euler characteristic (j2D is determined by the spectrum 
of Hk ■ Thus the direction of the logic in [16] is roughly the reverse of that in the present 
paper. Our derivation of Corollary [231 is simpler (as well as more general). 

6. Conclusions and additional remarks 

We have demonstrated that the "topological" term in the heat-kernel expansion for a 
Laplacian on a quantum graph does indeed have an index interpretation, if the Laplacian 
is of the scale-invariant class (i.e., the boundary conditions do not mix function values 
and derivatives, and hence the scattering matrix is independent of k). Such a Laplacian 
factors into two first-order operators, A and A*, defined on domains determined by those 
boundary conditions. 

One can calculate the index in three ways: (1) in the usual way, by subtracting the 
heat kernel of AA* from that of A* A (Proposition [16] and Theorem [201) : (2) by inspection 
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of just the heat kernel of A* A (Proposition [TBI and Corollary [TUj) ; (3) just by counting 
the number of Dirichlet-type conditions (Theorem [TH and Corollary [T5|) . 

A general index formula has been derived for an arbitrary (elliptic) differential 
operator on a quantum graph (Theorem [H]). 

Along the way, we have provided some properties of the Hamiltonian that are 
equivalent to the scale invariance of the vertex conditions (Theorem [HJ. 

Finally, we have determined the algebraic multiplicity of as a root of the secular 
equation of a generic scale- invariant graph Laplacian in a novel way (Corollary 1231. 

We now add a few final remarks concerning the results of the paper: 

• There is an elementary sense in which the integer V — E encountered in (T5]) is 
associated with an operator index. A graph as a purely combinatorial object (with 
no lengths assigned to the edges) is described in graph theory by the incidence 
matrix, whose rows are indexed by the vertices and its columns by the edges. Each 
matrix entry is equal to either 0, 1, or 2, depending on whether that edge is not or 
is attached to that vertex or forms a loop there. Then it is easy to see that (just 
because of the matrix's dimensions) the index of the incidence matrix is equal to 
E — V. 

• The restriction to "local" vertex conditions in Theorem [2] and elsewhere has very 
little content. Indeed, the structure of the graph enters the problem only through 
the vertex conditions, and one could define a vertex as a subset of edges that are 
related by such conditions. Alternatively, one could think that all vertices of a 
quantum graph have collapsed into a single one (creating a "rosette" consisting of 
one vertex and E cycles attached). Since all (whether previously local or nonlocal) 
vertex conditions refer to this single vertex, all conditions have become local. This 
is impossible only if we need to enforce some specific type of vertex conditions, e.g., 
the Kirchhoff ones; after the graph collapses to a rosette, the vertex conditions will 
generally no longer be of that type. 

• Kirchhoff conditions on A°(r) and anti-Kirchhoff conditions on A 1 (r) seem to be 
quite natural, whereas the interchanged conditions look rather unnatural. The 
situation for manifolds is different: there one can either impose Neumann conditions 
on forms of even degree and Dirichlet conditions on forms of odd degree, or do the 
reverse. The two choices correspond to two different cohomology theories for the 
manifold, "absolute" and "relative" [61 [7] . 

• Carlson [2] also constructs second-order self-adjoint operators on quantum graphs 
in terms of first-order operators, but his construction is rather different from ours. 
More pertinent are the remarks of Friedman and Tillich [5] and Exner and Post [3] 
that derivatives on quantum graphs should be treated as 1-forms (or vector fields). 

• An extension of the first-order formulation of the index theorem to operators H 
with a nontrivial Robin part, C v ^ 0, is not to be expected. As pointed out 
by V. Kostrykin, in general either H or its dual will have negative spectrum, 
whereas operators of form A* A and AA* must both be positive. This impossibility 
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of factorization is also contained in the statement (vii) of Theorem [HJ On the other 
hand, inserting "Robin" operators A„ into the definition of H does not change the 
t-independent terms of Tr Kh and its dual, which are the ingredients of the formula 
for the index of the second-order operator. 
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